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Abstract. We discuss the dynamics of an open two- mode Bose-Hubbard system 
subject to phase noise and particle dissipation. Starting from the full many-body 
dynamics described by a master equation the mean-field limit is derived resulting 
in an effective non-hermitian (discrete) Gross-Pitaevskii equation which has been 
introduced only phenomenologically up to now. The familiar mean-field phase space 
structure is substantially altered by the dissipation. Especially the character of the 
fixed points shows an abrupt transition from elliptic or hyperbolic to attractiv or 
repulsive, respectively. This reflects the metastable behaviour of the corresponding 
many-body system which surprisingly also leads to a significant increase of the purity 
of the condensate. A comparison of the mean-field approximation to simulations of the 
full master equation using the Monte Carlo wave function method shows an excellent 
agreement for wide parameter ranges. 
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The physics of ultracold atoms in optical lattices has made an enormous progress 
in the last decade, as it is an excellent model system for a variety of fields such as 
nonlinear dynamics or condensed matter physics [1,2]. Although this seems to be 
an inherent many-particle problem, the dynamics of the macroscopic wave function 
is remarkably well reproduced by the (discrete) Gross-Pitaevskii equation (GPE) if 
the systems undergoes a Bose- Einstein condensation [3]. Recently there has been an 
increased theoretical [4-8] as well as experimental [9] interest in the dynamics of these 
systems coupled to the environment. In particular, the effects of particle loss have been 
discussed in several theoretical articles resorting to an effective non-hermitian mean- 
field description introduced phenomenologically to analyse resonances, transport and 
localization effects [10-16]. 

In this communication we want to illuminate the origin and give a convincing 
motivation of this approach. Starting from a master equation describing the full many- 
body dynamics including phase noise and particle loss we derive a generalized, non- 
hermitian Gross-Pitaevskii equation. Due to the decay the structure of the resulting 
dynamics abruptly changes introducing repulsive and attractive fixed points. Unlike the 
dissipation-free case there are no longer oscillations around the fixed points, such that 
irreversible transitions between the former self-trapping fixed points are possible. This 
reflects the metastable behaviour of the open many-particle system and gives rise to a 
significant purification of the BEC by the dissipation which will be explained here. The 
validity of the presented approximation is tested by a comparison to full quantum many- 
body calculations, showing that the mean-field approximation provides an excellent 
description of the system. To integrate the effective description by a non-hermitian 
GPE into well-known concepts from the theory of open quantum systems, we discuss 
their relation to the quantum jump approach [17,18], demonstrating that they provide 
a well-suited tool to analyse the short- as well as the long-time behaviour of the open 
many-particle system. 

In particular, we consider the dynamics of ultracold atoms in an open double-well 
trap which is not only an extremely popular model system, but also has several recent 
experimental realisations [19-22]. The unitary part of the dynamics is described by the 
Bose-Hubbard type Hamiltonian 

H — — J (h\a,2 + a\aij + e\h\ + e 2 n 2 + ^ (a^a\ + a^afj , (1) 

where cL and a\ are the bosonic annihilation and creation operators in mode j and 
t 

hj = CLjCLj is the correponding number operator. We set H — 1, thus measuring all 
energies in frequency units. In order to analyze the dynamics in the Bloch-representation 
we transform the Hamiltonian using the collective operators 

L x = - (a\a 2 + , L v = \ {^\ a 2 ~ 4«i) > 

which form an angular momentum algebra su(2) with rotational quantum number 
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N/2 [6,23-25]. With these definitions the Hamiltonian ([I]) can be rewritten as 

H = -2JL X + 2eL z + UL 2 Z (3) 

with 2e = ei — e 2 up to terms which only depend on the total number of atoms. 
The macroscopic dynamics of the atomic cloud is well described within a mean- field 
approximation, only considering the expectation values of the angular momentum 
operators £j = (Lj) and the particle number n = {hi + h 2 ) [6,23,24]. 

Here we consider the dissipative extension of this system. The main source of 
decoherence in current experiments is phase noise due to elastic collision with atoms in 
the thermal cloud [4, 5] which effectively heats the system, leaving the particle number 
invariant. Only recently, methods to tame this source of decoherence were discussed 
in [30,31]. In this communication we focus on the effects of particle loss at constant 
rates j a j in the two wells j = 1,2. Such a loss is not only of fundamental interest but 
can be adapted in current experiments without greater difficulties by removing atoms 
with a focused resonant laser beam or by inducing a radio frequency transition to an 
untrapped internal state [26] . All parameters used here lie in realistic ranges for ongoing 
experiments [19,20]. 

The master equation description including both, phase noise and particle loss, is 
well established [27] and routinely used in the context of photon fields. Thus we consider 
the dynamics generated by the master equation 

p= - i[H, p] - ^2 (n 2 jP + P™j ~ Zrijphj) 

A j=l,2 

- 2 Yl ^ (f&P + ph Y a i ~~ 2 %P a j) • ( 4 ) 
i=i,2 

The evolution equations for the expectation values of the angular momentum operators 
(|2J) can be calculated starting from the Master equation via £j = tr(Ljp) with j = x, y, z. 
This yields the exact result 

£ x 2e£y 2U (£y£ z -\- Ay Z ^) 7~2 £ x ? 
£y = 2Jt z + 2e£ x + 2U(£ X £ Z + J X z) ^2 ^"y> 
£ z = 2J£y T-y £ z T-y f a n/2, 

h = - T^n - 2TfV a 4, (5) 

where we have defined the transversal T x _1 and longitudinal T 2 _1 damping rates by 

Tf 1 = (7ai+7a 2 )/2 and T^ 1 = 7 P + T-f 1 (6) 

and the asymmetry of the loss rates by /„ = (7 a2 - 7oi)/(7oi + 702) ■ I n the non- 
interacting case U = 0, these equations of motion ressemble the Bloch equations in 
nuclear magnetic resonance [28], except for the fact that the 'equilibrium' value of the 
population imbalance £ z is given by —f a n/2 and therefore depends on the decreasing 
expectation value of the total particle number n. 

The exact equations of motion (jSJ) still include the covariances = (LjLk + 
LkLj)/2 — (Lj)(Lk). The approximation of second order moments by products of 



Mean-field dynamics of a two-mode BEG subject to noise and dissipation 4 




Figure 1. Comparison of the mean- field approximation (thin blue line) with the 
full many-particle dynamics calculated with the MCWF method (thick red line) for 
,/ = 10s _1 , U = Is -1 , 7 P = 3s _1 and 7 a 2 = 5s _1 . The initial state was assumed to be 
a pure BEC (i.e. a product state) with s = (0.46,0,0.2) (a) and s = (0.14,0,0.48) (b), 
respectively, and n — 100 particles. 

expectation values, such that w yields the well-known mean-field description. 
This truncation is valid in the macroscopic limit, since the covarinces vanish as \ jn if 
the many-particle quantum state is close to a pure BEC. Here and in the following we 
depict the rescaled variables Sj = £j/n, thus renormalizing to separate the decay of the 
particle number n from the internal dynamics. 

The benefit of the mean-field approximation is illustrated in figure [I], where it 
is compared to the full many-particle quantum dynamics calculated with the Monte 
Carlo wave function (MCWF) method [17, 18]. The trajectory in figure [1] (a) was 
launched at s = (0.46, 0, 0.2) with a moderate population imbalance, thus performing 
Josephson oscillations [19]. The amplitude is damped because of the phase noise, 
while the oscillation period increases as the effective macroscopic interaction strength 
g(t) := Un{t) decreases. The decay of the particle number n(t) is also strongly 
modulated by the oscillations of the population imbalance. The trajectory in figured] (b) 
was launched at s = (0.14,0,0.48) in the self-trapping region. The residual oscillations 
are rapidly damped out and the system relaxes to a quasi-steady state on the shown time 
scale. The particle number decreases slowly and non-exponentially, since the condensate 
is mostly localized in the non-decaying potential well, cf. also [15]. All these features 
of the dynamics are well reproduced by the mean-field description, and the decay of 
the particle number is accurately predicted. Strong deviations are only expected in the 
vicinity of unstable fixed points of the mean-field dynamics, which can be nearly cured 
within the framework of phase space distributions [25]. 

In order to explore the genuine effects of particle loss, phase noise is neglected 
(Tp = 0) in the following. In this case the dynamics can be further simplified and one 
can easily show that the particle number coincides with the magnitude of the Bloch 
vector -)- £2 _|_ £2 _ 71/2, which we can use to reformulate the mean- field dynamics 
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Figure 2. Mean- field dynamics for the non-interacting case 5 = (upper row) and 
for a fixed interaction strength g = 4s _1 (lower row) in dependence of the decay rate 
(7 = for (a) and (d), 7 = 1.9 s _1 for (b), 7 = 2.1 s -1 for (c), 7 = Is -1 for (e) and 
7 = 4s _1 for (f)) - for all figures holds J = 1 s _1 and e = 0. 

by an effective non-hermit ian Gross-Pitaevskii equation 

• d_ ( 1>i\ = ( h + U\^\ 2 -J (7) 
dt\1> 2 ) { -J i 2 + U\ij 2 \ 2 J \ij 2 J {) 

with complex on-site energies ij = ej — ijaj/2. The equivalence to the Bloch vector 
description is established via the identification 

«. = ^| a - W 2 ) (8) 

and n = IV'iI 2 + |V ; 2| 2 - I n this effective description a loss of particles is represented by a 
loss of normalization. 

We now consider the dynamics for a fixed value of the macroscopic interaction 
strength g = Un = const, in the special case j a i =: 7 and 7 a 2 = 0. Even though the 
restriction to a fixed interaction constant seems to be artificial, it reveals the effects of the 
particle loss on the structure of the mean-field phase space and especially the character of 
the fixed points most clearly. Moreover, the dynamics under a fixed interaction constant 
correspond to the periods of constant particle number between two loss processes in 
the quantum jumps picture [17,18]. Therefore this treatment provides a well-suited 
description of the short- as well as the long-time behaviour. Note that the more general 
case j a i 7^ and •y a 2 7^ does not lead to a fundamentally different dynamical behaviour 
since only the difference of the decay rates influences the internal dynamics. However, 
the expectation value of the particle number n and thereby also the effective interaction 
strength g(t) decrease faster. 
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The resulting dynamics of the Bloch vector is illustrated in figure [2j The upper 
row (a-c) shows the phase space for the linear case, U = 0, where the mean-field 
approximation is exact. Without loss one recovers the famous Josephson oscillations 
(a). An analysis of the fixed points for the dissipative dynamics shows the emergence 
of two regimes depending on the amplitude of the loss rate. For weak losses, |7| < 2 J, 
the fixed points are given by 

/ ±[i-(£) 2 ]*\ 



SL 
4J 



V 







(9) 



/ 



While the fixed points remain elliptic and the population is still equally distributed, the 
fixed points are no longer symmetric, since the relative phase between them decreases 
(b). This behaviour can be qualitatively understood within the analogy to Josephson 
junctions: The weak decay induces an assymetry between the wells leading to to a 
continuous particle stream to the first well. At the fixed points this effect is compensated 
by the Josephson current Ij oc Js y requiring s y 7^ 0. 

For stronger decay rates, I7I > 2J, the two fixed points are given by 

/ \ 

(10) 
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Above the critical value I7I = 2J the character of the two fixed points changes abruptly 
from elliptic into an attractive and a repulsive one as shown in figure [2] (c). The maximal 
Josephson current is no longer sufficient to compensate the current induced by the decay 
leading to a population excess in the non-decaying site. This explains the population 
imbalance in the fixed points which increases with ascending decay rates. 

In the strongly interacting case without dissipation one observes the splitting 
of one of the elliptic fixed points into two novel elliptic and one hyperbolic fixed 
point - this is the famous self-trapping effect [19,23,24]. The critical interaction 
strength for the occurence of this bifurcation is lowered in the presence of dissipation to 



= U 2 n 2 > 4J 2 
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In the subcritical regime for 7 < 2 J and Un < 4 J 2 — 7 2 , we find 



oscillations around the same fixed points s± as in the non-interacting, but dissipative 
case ([9]). However, these are now distorted (not shown in the figure). In the overcritical 
regime g = Un > \J '4 J 2 + 7 2 and for a weak decay 7 < 2 J one rediscovers a generalized 
self-trapping effect. As a result of the dissipative process, one elliptic fixed point now 
bifurcates into an attractive and a repulsive fixed point (in contrast to the two elliptic 
ones) and one hyperbolic one (cf. figure [2] (e)). The novel fixed points are located at: 
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Figure 3. Comparison of the many-particle dynamics (thick red line) to the mean- 
field approximation (thin blues line) for an initially pure BEC. with s = (—0.5,0,0) 
and n(0) = 200 particles: Dynamics of the Bloch vector s (a), evolution of the purity 
of the BEC (b) and evolution of the population imbalance s z (c). The system initially 
relaxes to a non-linear quasi-steady state with a purity of almost one which is then 
lost as n(t) decreases. Parameters are chosen as J = Is -1 , Un(0) = 10s _1 , T\ = Is 
and f a = 1. 

For stronger decay rates, 7 > 2 J, the hyperbolic and the elliptic fixed point s± ([9]) 
meet and annihilate themselves as illustrated in figure [2] (f). Their disappearence is 
accompanied by the complete disintegration of periodic orbits. 

Let us finally discuss the implication of this phase space structure. We especially 
focus on the emergence of the attractive fixed point since it is stable and therefore 
strongly influences the many-body quantum dynamics. Figure [3] shows the dynamics of 
the rescaled Bloch vector s comparing results of a MCWF simulation (solid red line) 
to the mean- field approximation (thin blue line). The given parameters correspond to 
the situations illustrated in figure [2] (b) or (e), respectively, depending on the value of 
the macroscopic interaction strength g(t) = Un(t). The Bloch vector first relaxes to 
the attractive fixed point illustrated in figure [2] (e). The contraction of the mean- field 
trajectories to the attractive fixed point manifests itself by a convergence towards a pure 
BEC, which is the state of tightest localization in phase space [25]. This is illustrated 
in figure [3] (b) where we have plottet the purity V : = 2tr(p^ ed ) — 1 of the reduced single 
particle density matrix p re d, V = 1 indicating a pure BEC [6]. However, the attractive 
fixed point is lost as g(t) = Un(t) decreases, and thus the Bloch vector departs again. 
This behaviour is very well prediced by the mean-field approximation already for the 
modest atom number in the simulation. The mean-field trajectory then convergences 
to the limit cycle shown in figure [2] (b). However, as the atoms are so rapidly lost 
nearly no particles remain to follow the limit cycle predicted by mean-field theory. This 
transition effect between different fixed points is closely related to the quantum state 
diffusion in and out of a metastable state, which can be observed in optical bistability 
(see, e.g., [32]). Note however, that the system considered here irretrievably departs 
from the metastable self-trapping state because the fixed point is lost as n(t) decreases. 

In summary, we have derived a mean-field approximation for a dissipative two- 
mode BEC, starting from the full many-body dynamics described by a master 
equation including phase noise and particle losses. This treatment puts the so 
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far phenomenological description of open systems via non-hermitian Gross-Pitaevskii 
equations on a firm footing and paves the way for a variety of future applications in 
particular because the extension to an arbitrary number of modes is straightforward. 

An analysis of the resulting equations for a fixed interaction contant g shows that 
not only the critical value for the self-trapping bifucation is lowered but also the character 
of the fixed points abruptly changes to attractive and repulsive, such that one of them 
becomes unstable. Taking into account the decline of the interaction constant g(t) due 
to the particle losses, the system initially converges to the attractive fixed point, but 
then suddenly jumps to a Josephson oscillation as soon as g(t) falls below the critical 
value. This effect is understood as a manifestation of the metastable behaviour of the 
many-particle system and leads to a significant increase of the purity of the quantum 
state compared to the dissipation-free case. 

The comparison to numerical results for the many-particle system obtained via 
the Monte Carlo Wave function method shows that the approach presented here is an 
excellent approximation to the full many-body dynamics already for a modest initial 
number of atoms. Thus it provides an excellent basis for a further analysis of the 
interplay between dissipation and interaction [31]. Likewise the embedding into the 
more general concept of mean-field description [33] will be subject to future research. 

Acknowledgments 

This work has been supported by the German Research Foundation (DFG) through 
the research fellowship programm (grant number WI 3415/1) and the Heidelberg 
Graduate School of Fundamental Physics (grant number GSC 129/1), as well as the 
Studienstiftung des deutschen Volkes. We thank H J Korsch, E M Graefe and A Niederle 
for stimulating discussions. 

References 

[1] Morsch O and Oberthaler M 2006 Rev. Mod. Phys. 78 179 

[2] Madronero J et al. 2006 Adv. At. Mot Opt. Phys. 53 33 (ed M Scully and G Rempe) 

[3] Pitaevskii L and Stringari S 2003 Bose-Einstein Condensation (Oxford University Press). 

[4] Anglin J R 1997 Phys. Rev. Lett. 79 6 

[5] Ruostekoski J and Walls D F 1998 Phys. Rev. A 58 R50 

[6] Vardi A and Anglin J R 2001 Phys. Rev. Lett. 86 568; Anglin J R and Vardi A 2001 Phys. Rev. 
A 64 013605 

[7] Ponomarev A V, Madronero J, Kolovsky A R and Buchlcitncr A 2006 Phys. Rev. Lett. 96 050404 
[8] Wang W, Fu L B and Yi X X 2007 Phys. Rev. A 75 045601 

[9] Syassen N, Bauer D M, Lettner M, Volz T, Dietze D, Garcia-Ripoll D D, Cirac J I, Rempe G and 
Diirr S 2008 Science 320 1329 
[10] Moiseyev N and Cederbaum L S 2005 Phys. Rev. A 72 033605 
[11] Schlagheck P and Paul T 2006 Phys. Rev. A 73 023619 

[12] Witthaut D, Graefe E M, Wimberger S and Korsch H J 2006 Phys. Rev. A 75 013617 
[13] Livi R, Franzosi R and Oppo G L 2006 Phys. Rev. Lett. 97 060401 
[14] Hitler M, Kottos T and Ossipov A 2006 Phys. Rev. A 73 063625 



Mean-field dynamics of a two-mode BEC subject to noise and dissipation 



9 



[15] Schlagheck P and Wimbergcr S 2007 Appl. Phys. B 86 385 

[16] Graefe E M, Korsch H J and Niederle A 2008 Preprint larXiv:0807.1777l 

[17] Dalibard J, Castin Y and M0lmer K 1992 Phys. Rev. Lett. 68 580 

[18] Carmichael H J 1993 An Open Systems Approach to Quantum Optics (Springer, Berlin) 
[19] Albiez M, Gati R, Foiling J, Hunsmann S, Cristiani M and Oberthaler M K 2005 Phys. Rev. Lett. 
95 010402 

[20] Gati R, Hemmerling B, Foiling F, Albiez M and Oberthaler M K 2006 Phys. Rev. Lett. 96 130404; 
Gati R, Esteve J, Hemmerling B, Ottenstein T B, Appmeier J, Weller A and Oberthaler M K 
2006 New J. Phys. 8 189 

[21] Schumm T, Hofferberth S, Andersson L M, Wildermuth S, Groth S, Bar- Joseph I, Schmiedmayer 

J and Kriiger P 2005 Nature Physics 1 57 
[22] Foiling S, Trotzky S, Cheinet P, Feld M, Saers R, Widera A, Miiller T and Bloch I 2007 Nature 

448 1029 

[23] Milburn G J, Corney J, Wright E M and Walls D F 1997 Phys. Rev. A 55 4318 

[24] Smerzi A, Fantoni S, Giovanazzi S and Shenoy S R 1997 Phys. Rev. Lett. 79 4950 

[25] Trimborn F, Witthaut D and Korsch H J 2008 Phys. Rev. A 77 043631; Trimborn F, Witthaut D 

and Korsch H J 2008 Preprint arXiv:0802.1142l 
[26] Bloch I, Hansen T W and Esslinger T 1999 Phys. Rev. Lett. 82 3008 
[27] Gardiner C W and Zoller P 2004 Quantum Noise (Springer Series in Synergetics, Berlin) 
[28] Bloch F 1946 Phys. Rev. 70 460 

[29] Tikhonenkov I, Anglin J R and Vardi A 2007 Phys. Rev. A 75 013613 

[30] Khodorkovsky Y, Kurizki G and Vardi A 2008 Phys. Rev. Lett. 100 220403 

[31] Witthaut D, Trimborn F and Wimberger S 2008 in preparation 

[32] Rigo M, Alber G, Mota-Furtado F and O'Mahony P F 1997 Phys. Rev. A 55 1665 

[33] Duffield N G and Werner R F 1992 Rev. Math. Phys. 4 383 



